TUM-HEP 215/95 
IC/95/29 
hep-ph/95"0323^ 

February 1995 



Axino Mass in Supergravity Models 



E. J. Chun* 

International Center for Theoretical Physics 
P. O. Box 586, 34100 Trieste, Italy 

and 

A. Lukas^ 

Physik Department 
Technische Universitat Miinchen 
D-85747 Garching, Germany 
and 

Max-Planck-Institut fiir Physik 
Werner-Heisenberg-Institut 
P. O. Box 40 12 12, Munich, Germany 



Abstract 

We analyze the mass of the axino, the fermionic superpartner of the axion, in general 
supergravity models incorporating a Peccei-Quinn-symmetry and determine the cosmo- 
logical constraints on this mass. In particular, we derive a simple criterion to identify 
models with an LSP-axino which has a mass of 0{m^i^l fpQ) = O(keV) and can serve 
as a candidate for (warm) dark matter. We point out that such models have very special 
properties and in addition, the small axino mass has to be protected against radiative 
corrections by demanding small couplings in the Peccei-Quinn-sector. Generically, we 
find an axino mass of order 772-3/2 • Such masses are constrained by the requirement of 
an axino decay which occurs before the decoupling of the ordinary LSP. Especially, for 
a large Peccei-Quinn-scale fpq > 10^^ GeV this constraint might be difficult to fulfill. 
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The implications of axions have been examined extensively since their existence was sug- 
gested by an attractive mechanism for resolving the strong CP problem ^ . Even though 
axions are very weakly interacting, their astrophysical and cosmological effects are strong 
enough to narrow down the window of the Peccei-Quinn-scale fpq to 10^° GeV < fpg < 
10^^ GeV P). On the same footing the axino as the supersymmetric partner of the axion can 
play an important role in astrophysics and cosmology . An interesting feature is that axinos 
may receive a mass of order keV which would render them a good candidate for warm dark 
matter. If axinos are heavier than a few keV they have to decay fast enough not to upset 
any standard prediction of big-bang cosmology. Given the weakness of their interactions, a 
constraint on their lifetime put a rather severe limit on the lower bound of their mass. There- 
fore it is very important to know the axino mass in discussing the cosmological implications 
of supersymmetric axion models. 

In global supersymmetry (SUSY) the calculation of the axino mass was performed in 
refs. [0, §]. In this paper, we will provide the computations in models with local supersym- 
metry (supergravity). Some partial results have been obtained in refs. |P, |T0|- In the case 
of spontaneously-broken global SUSY the axino mass is of the order Tn\i2/ fpq ~ keV where 
m3/2(~ ITeV) is taken to be the global SUSY-breaking scale |^. On the contrary, in the 
context of supergravity, the axino mass can be of order m^j2 as first noticed in ref. 0. Soon 
after this it was realized that the axino mass is truly model-dependent and the global SUSY 
value rn\i2/ fpQ niay be obtained in supergravity models as well |T^. We will extend those 
results in a generic treatment of supergravity models also including radiative corrections. 

The prime motivation for supergravity is well-known. Realistic supersymmetric general- 
izations of the standard model are based on local SUSY spontaneously broken in a so-called 
hidden sector at a mass scale of order Ms ~ 10^^ GeV [0]. The induced SUSY breaking 
scale in the observable sector is determined by a value of the order of the gravitino mass 
77^3/2 ~ Mg/Mp where Mp is the Planck scale. Axionic extensions of the minimal super- 
symmetric standard model (MSSM) inevitably incorporate an extra sector which provides 
spontaneous breaking of the Peccei-Quinn [/(l)-symmetry at the scale fpq. This sector 
(PQ-sector) is considered as a part of the observable sector. In the framework of effective 
supergravity theories with a Lagrangian composed out of a global SUSY part and soft terms 
the hidden sector dependences are encoded in the soft terms. We will rely mostly on this 
effective approach as it makes the calculations tractable. 

A color anomaly in the PQ-sector can be introduced in two ways. The fields S in this 
sector can be coupled to the standard Higgs doublets Hi, H2 of the MSSM like gSHiH2 
(DFSZ-axion) ^ or to new heavy quarks Qi, Q2 like gSQiQ2 (KSVZ-axion) [Q. Our 
analysis of axino mass will be concerned with the tree level result in the effective theory which 
is obtained after breaking the PQ-symmetry and SUSY and should therefore not depend on 
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which implementation of the axion is chosen. We will find that this mass crucially depends 
on the structure of the PQ-sector as well as on the hidden sector. It should, however, be 
mentioned that in addition the above couplings and the corresponding trilinear soft terms 
lead to a one-loop radiative mass of order m^ ioop — 3g'^Am^/2/^Q'^'^ is the trilinear soft 
coupling) in], p. For the DFSZ-axion this contribution is clearly negligible since the coupling 



g has to be quite small {g ~ 10^'') in order to generate reasonable Higgs-masses. In the 
KSVZ-implementation there is no such restriction on g and the above one-loop contribution 
can modify the result for the axino mass to be obtained below. 

An interesting observation is that the axino mass depends on whether the PQ-sector 
admits additional accidental zero modes on the global SUSY level. As a general statement we 
find that the axino mass cannot be bigger than 0(7713/2) if the axion mode is the only global 
zero mode in the PQ-sector. In order to see this, it is useful to analyze the full supergravity 
Lagrangian. Later, we will confirm this result by using the effective Lagrangian approach. 

The PQ-sector consists of an arbitrary number of singlets S"^ with charge qa under the 
PQ-symmetry and provides its spontaneous breaking at the scale fpg = {Y,a lll^^l'^Y^'^^ where 
_ (^g<iy The axion multiplet $ is given by $ = J2a QaV'^S"' / fpq. Its component field content 
reads $ ~ {s + ia, a) with the axion a, the saxion s and the axino a. Supersymmetry is broken 
mainly by the hidden sector with singlet fields Z*. To simplify the argument we assume 
minimal kinetic term for the PQ fields as well as for the hidden sector fields. The scalar 
potential reads 



V = M], exp[G/M|,] ( ^ G,G* + ^ GaGl - 3M| 



where G = K + Mphi\W/Mp\^ and W = W{S'') + W{Z'). Spontaneous supersymmme- 
try breaking implies (Gi) ~ Mp for some i and generation of a gravitino mass, m3/2 — 
Mpexp[G/2Mp]. We also have to admit the possibility that (Go) ~ Mp for some a. We now 
look at the minimization condition (Vi) = to estimate the axino mass. Vanishing of the 
cosmological constant (V) = is assumed. Then we have 

(K) = {GabGl + Ga^: + Gl) . (2) 

If we take a massive mode S"" in the PQ-sector, {Gab) is dominated by 

{Gab)^{Ml^)^I^6ab. (3) 

W 1713/2 

since (W) ~ M|,m3/2. In addition, we can estimate the maximal order of magnitude of (Gai) 

- i-Ml'^} ^ ^ (4) 
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since {Wi) ~ Mpm^/2 maximally and {Wa) ~ fpQ- Therefore the condition (K) — 
gives {Ga) ~ fpQ- In the axion direction — JZaQaV'^Ga/ fpQ, we also have (G^) — 
{Y.aQaV"'Ka/ fpq) — Y^a Qi\'^°'\'^ / fpQ ~ fpQ- Sincc the absence of accidental zero modes is 
required in the PQ-sector one finds that (Ga) ~ fpQ for all the fields S"'. 

It is now straightforward to estimate the maximal order of the axino mass by using the 
fermion mass matrix in supergravity models, Mij— Mpexp{G/2Mp)[Gab + GaGb/Mp]. Along 
the axino direction, 

From the U{1) invariance of G, we get J2aQaVa{Gab) = Qbivb — {Gb)) ~ fpq- Hence the axino 
mass is maximally of order 1713/2- 

The above order-of-magnitude estimation is indeed insensitive to the specific forms of 
the kinetic term or the supcrpotential as long as the hidden sector fields have only non- 
renormalizable couplings to the observable sector. This happens because higher power terms in 
KorW are naturally suppressed by powers of Mp which renders their contribution negligible. 
Therefore we conclude that the axino mass in general supergravity models is at most of the 
order 777.3/2 if no other zero mode than the axion is present in the PQ-sector. 

On the other hand if there are extra zero modes we are not able to constrain further the 
order of (Go) for a zero mode direction a so that the above argumentation breaks down. In 
fact, axino masses ^ 1713/2 are possible in those models as we will see below. 

We will now calculate the actual value of the axino mass relying on the effective super- 
gravity Lagrangian with arbitrary soft terms. This allows to consider axino masses down to 
O {ml/2/ fpq) where the next to leading order in the 1/Mp-expansion of supergravity becomes 
important. 

The supcrpotential W of the PQ-sector is expanded as 

W = /a5« + ^fabS'^S'' + ^fabcS'^S'S'^ (6) 

and a departure from standard soft terms is encoded in 

N = + ^dabfabS^^s' + ^dabJabcS^'s's' . (7) 

Then the scalar potential reads 



V = d'WdbW + mfsbs" + mg/s [s'dbW + {A-?,)W + N + h.c.J . (8) 

To minimize this potential we apply the following strategy : The VEVs v"" are split into a 
global SUSY value u"- with daW{u) ~ and corrections due to the soft terms, v"" — u°'+w"'. 
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Expanding daV around the global minimum (u") results in 

idaV){v) = ^fabeP'''w,WdW'+^MatP'''w,Wd+M'j 

+ ^{A + dabc)fabcW'w' + 771^2 [((A - 1)M„5 + dabfab)w' (9) 



+ (1 + dabc)fabcU''w'\ + mfwb + M^bM^'w, + la 

with the global mass matrix 

Mab = dadbW{u) (10) 

and the definitions 

la = my^Ja + ml'ub , Ja = Mabu' + {daN){u) . (11) 

Care should be taken on the choice of (m"). With any global minimum also a rotation 
exp{qaz)u"' with z = x + iy under the complexified Upq{1) leads to such a minimum. Despite 
the ^/-dependent part of this symmetry which is clearly present in the full theory the x- 
dependent part (present because of the holomorphy of the superpotential) is broken by the 
soft terms. Therefore an appropriate fixing for the x-dependent part of the symmetry should 
be applied such that the global minimum comes close to the local values (w") resulting in small 
expansion coefficients (w°). Such a fixing is provided by the condition Wa = J2a QaU^'w"' / fpQ = 
implying that the correction in the global axion direction (denoted by an index a) vanishes. 
As can be expected the axino mass is expressible in terms of the corrections (w") : 

iM~a)a = -^qt{dbW){v) 



^ q\ [Ubcw' + ^fbcdW^w'' + ■ ■ ■) . (12) 



These corrections have to be determined from eq. (|^). Let us work in a basis with diagonal 
global mass matrix Mab = MaSab- We denote massive modes with indices ■ ■ ■ and possible 
additional zero modes with indices /5, 7, ■ ■ ■. Then an important observation is that the cor- 
rections Wi are basically determined by the linear term |Mjpw* and Jj in eq. (y). In zero mode 
directions the situation might be more complicated since |M^pw^ can be small compared to 
other terms in eq. (^. Taking this into account we conclude that the expression 

^ "^^1^^"' ~ ife^"^^^'''' " 2^^'^"^^'^'''^" + 0{mljJfpQ) (13) 

gives the correct order of the axino mass. For the moment we leave the values of unspecified. 
Instead we concentrate on the first term in eq. (|T^) and spht the expression into its Planck 
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(or GUT)-scale value J^^-* and corrections J^^-* arising from renormalization down to /pg. The 
generic value of the corrections can be roughly estimated as 

4" = KkJlq (14) 
K^3.^nm (15) 



with appropriate combinations Aa, of the superpotential couplings. According to eq. (|TT]) 
J^°^ is naively of the order fpQ. Therefore a necessary condition for the axino mass to he much 
smaller than m3/2 is that j'f'^ = 0. 

In general this condition implies a relation between the structure of the superpotential and 
the soft terms. It is instructive to analyze this relation for a certain subclass of models, namely 
those with independent soft coupling A, B, C for the trilinear, bilinear and linear terms in the 
superpotential, respectively. A computation leads to 

= (5 - C)Mabv!' + ^(A -2B + C)fabcv!'u' . (16) 

Depending on the properties of the superpotential (and assuming that at least one coupling 
fa is nonzero to force the symmetry breaking) several cases can be distinguished : 

• fabU^-i fabcU^u'^-i MabU^ 7^ : Then J^°^ = if and only ii A = B = C . No special property 
of the superpotential is required. 

• fabU^ = 0, Mu 7^ : Then J^'^^ = if and only li A = C. A simple superpotential which 
fulfills this requirement is e. g. W = X{SS' — fi'^)Y since all fab = 0. 

• fabcU^u" = 0, Mm ^ : Then ) = if and only if B = C. 

• Mu = 0, fabv!' i- : Then jf^ = if and only if A - 25 + C = 0. As the only ones 
these models allow for the full standard pattern B = A — \^C = A — 2. They possess, 
however, at least one additional zero mode Yl,a'^"'^°'- An example is provided by the 
superpotential W = X{SS' - Z^)Y - X'{Z - jj^f 0. 

Observe that in particular for A = B = C the expression 7^°-' vanishes in any model. 

If no additional zero mode is present the axino mass is already completely determined by 
the first term in eq. (pIS]). On tree level this means 



JW = - mf^=0{ml/JfpQ) 



(17) 



We have therefore found a simple criterion to decide about the magnitude of the axino mass 
which for soft terms specified by the couphngs A,B,C singles out the particularly simple 
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patterns listed above. We remark that the small axino mass in no-scale models observed in 
ref. can also be understood in terms of our analysis since in those models A = B = C = 0. 
A full supergravity computation of the axino mass in no-scale models shows that their tree- 
level mass even vanishes. Therefore the first line of eq. (|17[) has to be understood as a generic 
result which in certain special cases might be too large. In the first case of a light axino mass 
radiative corrections to the potential parameters can become important. Using eq. (|14D this 
leads to a contributions of 

rria ~ fcim3/2 . (18) 

To keep the order rria = Oim^j^/ fpq) an upper bound on the couplings has to be required. 
For — 10^ GeV and fpq ~ 10^^ GeV this implies Aj ~ 10^^. A systematic way to avoid 
such small couplings is to consider no-scale models. Since gaugino masses are the only source 
of SUSY-breaking in those models the standard model singlet fields in the PQ-sector will not 
receive any radiative soft terms. 

We see that without additional zero modes a complete answer can be given. In particular 
we recover the result rria ~ 0(m3/2). 

If additional zero modes are present the situation becomes more complicated since the 
corrections wp in the zero mode directions can become large. In addition we have to consider 
that the zero entry Mp of the mass matrix receives a radiative contribution Mp = X'kpfpq 
with kp in analogy to eq. ( |T5| ) unless it is protected by an additional continuous or discrete 
symmetry. Let us discuss some relevant cases. First we discuss a model with jf^ = 0, e. g. a 
model of the last type in the above list for A, B, C-type soft terms. Then a tree level axino 
mass Oirni^i^/ fpq) is not guaranteed as opposed to the case without additional zero modes : 
If the terms in eq. (|^) linear in wp vanish (which e. g. occurs for /^j^ = 0, fjs^s = or A = 5) a 
value wp = 0(1^^^) results which causes an axino mass given by mf^ = O ((7/13/2 //pq)"'^^^"^3/2) • 
In any case a small axino mass <C /«3/2 has to be stabilized against radiative corrections. For 
superpotential couplings A = 0(1) the axino mass is shifted to rria = 0(7713/2). 

Now we assume that J^°^ 7^ 0. Then for values A = 0(1) the linear term |AfgpW/3 in eq. (^ 
will dominate and the axino mass is given by 

ma = 0{m3/2/kf3) . (19) 

1 /3 

If, on the other hand, the couplings A are very small we can have wp = 0{lp ) leading to an 
axino mass rria = 0{{m\i2fpqY^'^). We see that axino masses ^ rn^/2 are indeed possible. An 
example featuring all these aspects is the superpotential in the fourth entry of the above list 
for A, B, O-type soft terms. 

Now we turn to a discussion of the cosmological constraints on the masses of the axino and 
the saxion. We begin by noticing the fact that self-couplings among the axion supermultiplet 
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arise after integrating out the heavy fields in the PQ-sector 



-terms 



'1 + ^y^s) (^^df^ad^a + ^d^sd^s + ia'^^.d^a^ - -^df^aa-f^^a H , (20) 



where x = qfvf / fpQ . In some cases, in particular in a model with superpotential W = 
X{SS' — ^'^)Y and universal scalar soft masses, x is zero at the Planck scale and receives a 
contribution x ~ ln(Mp//pQ)/647r^ when the RG-improved potential at the PQ-scale is 



considered |13|. Generically, however, x is of order 1. In this case the self-coupling becomes 
important since a saxion can decay into two axions faster than e.g., into two gluons. Decay- 
produced axions do not heat the universe. Therefore the cosmological effect of saxion decay is 
different from what has been investigated assuming vanishing x [p!4| , |15|. If x is of order 1, a 
stronger bound on the saxion mass can be expected since the decay-produced axions are not 
thermalized but red-shifted away. The standard nucleosynthesis constrains the energy density 
of the universe due to this red-shifted axions to be less than what is contributed by one species 
of neutrinos at the time of nucleosynthesis. This gives the constraint msYgg^u < where 
Td = 0.55g^^^^\/TWp is the decay temperature of the saxion and F = x^m^/87r/pQ its decay 
rate. The relativistic degrees of freedom at are counted by g^^. We get 

This bound on the saxion mass which receives a contribution 0(m3/2) from scalar soft masses 
might be difficult to fulfill for large values of fpg. 

Cosmological implications of axinos were first discussed in ref. p assuming unbroken R- 
parity. The axino mass can be constrained in two ways. First, the axino can be the lightest 
supersymmetric particle (LSP). Then it should be lighter than a few keV in order not to 
overdose the universe. Otherwise, the axino decays into at least one LSP composed out of the 
neutralinos in the MSSM. In this case, the decay-produced neutralinos tend to overdominate 
the evergy density of the universe. To avoid this the axino should be heavy enough to decay 
before the neutralinos decouple. Considering the axino decay into photino plus photon, it was 
obtained that the axino mass should be bigger than a few TeV p. Then, the axino decay into 
top quark and scalar top can be allowed. From this decay channel one finds a less restrictive 
bound 

Here Xf is the PQ-charge of the top quark. We see that a wide range of axino masses between 
C»(keV) and 0(10^ GeV) is excluded. 
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Let us now analyze how the above constraints modify if an inflationary expansion is taken 
into account. The decoupling temperature of axino or saxion is around the range of the 
reheating temperature Tr ~ 10^° GeV which is the maximally allowed value to cure the 



gravitino problem in supergravity models [16]. If decoupling of the axino occurs before inflation 
the primordial axino relics are diluted away. The above consideration, then, has to be applied 
to the regenerated population of axinos. We recall that the axino decoupling is determined 
by its interactions with gluinos, quarks and anti-quarks via gluon exchange The axino 
decouples at the temperature 



The regenerated number density per entropy is given by 

y^2xlO-M^ (—-^\. (24) 



fpQ ) VlQiOGeV 

Depending on the range of the axino mass we can distinguish three cases as follows. 

First, for a stable axino, the constraint from overclosure gives the following loose bound 
on the axino mass in terms of Tr : 

Second, an axino with a mass satisfying the lower bound in eq. (22) is still allowed by cos- 



mology. Finally, for an unstable axino with mass between the estimations in eq. (|25D and 
in eq. (|2^), one gets a bound on the reheating temperature by replacing the axino mass in 
eq. (^Sf ) by the mass of the usual LSP since the decay-products of the axino contain at least 
one LSP: 

r«<2xlO»Gevf-#^Vfi^) . (26) 

VlO^GeVy \ m^o j ^ ^ 

This represents a quite stringent bound on Tq. 

In this letter we have analyzed the axino mass in general supergravity models and the 
cosmological constraints on such models. We have distinguished models with and without 
additional zero modes in the PQ-sector. For the latter we found the axino mass to be at 
most of 0(1713/2) ■ In the context of an effective approach encoding supersymmetry breaking 
in soft terms we were able to derive a simple necessary criterion for a small tree level axino 
mass <^ ^3/2 given in terms of superpotential and soft term properties. For uniform trilinear, 
bilinear and linear soft couplings A,B,C the criterion is always fulfilled for A = B = C whereas 
for A = B^C or Ay^B = C additional properties of the superpotential had to be required. 
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If the global vacuum {u) represents an additional zero mode of the globally supersymmetric 
theory, i. e. Mu = with the global mass matrix M, the relation A — 2B + C = which 
admits the standard pattern B = A — 1,C = A — 2 is sufficient for the criterion to hold. 

We showed that in models without additional zero modes our criterion is sufficient, i. e. it 
guarantees a tree level axino mass of at most mf^ = 0{m\i2/ fpo)- In the presence of other 
zero modes it serves as a good indication for such a small mass but additional conditions (like 
e. g. A ^ B for models with an A, B, C-pattern) are required to have mf^ = 0{my2/ J'pq)- 

From the cosmological point of view axino masses 0(171^2/ Ipq) = O(keV) are very attrac- 
tive. In this case the axino is the LSP and can contribute a relevant part of the mass in the 
universe as (warm) dark matter. Though models with such an axino mass can be constructed 
as we have seen they correspond to very special points in the space spanned by the superpo- 
tential and soft term parameters. Moreover, such small masses are not stable under radiative 
corrections arising from renormalization effects between Mp and the PQ-scale fpg. Taking 
these effects into account the axino mass will be generically given by m^^'' = 0{km^/2) with 
k = A^ln(Mp//pQ)/327r^ and a typical superpotential coupling A. The one-loop contribution 
from the characteristic coupling of the axino to the Higgs-fields or heavy quarks will be given 
by rria^ioop — 3g'^Am^/2/lQn'^ which is only relevant in the KVSZ-implementation of the axion. 
As the LSP the axino has to be lighter than a few keV which in turn puts a severe limit on 
the couplings A (and g in the KVSZ-case), typically A ~ 10^"^. If the decoupling temperature 
of the axino is larger than the reheating temperature of inflation the overclosure bound on 
the regenerated axino population is weakened resulting in a somewhat weaker bound on A, 
typically A ~ 10^'^. In any case we conclude that a cosmological relevant LSP-axino - though 
possible in principle - is not very likely to occur : Special models are needed and in addition 
small couplings have to be chosen in order to avoid a conflict with the overclosure bound. 

At this point it should be mentioned that no-scale supergravity models can provide a 
naturally light axino [Q. Those models are characterized by a special pattern of the soft 
terms : The only non-vanishing soft terms at the Planck scale are gaugino masses and therefore 
A = B = C = Oa.t tree level. Applying the above statements a small axino mass results in this 
case. In fact, a full supergravity calculation shows that the mass vanishes on tree level. Other 
soft terms for the gauge non-singlets can be generated due to renormalization group effects 
below the Planck scale in those models. Since the PQ-sector consists of singlets their soft 
terms are not affected by renormalization effects and the axino (saxion) remains massless 0. 
Non-vanishing masses can be however generated via the one loop contribution rria^ioop- In 
the DFSZ-implementation they are so small that cosmological effects of the axino and the 
saxino are negligible. This is clearly different in the KVSZ-case. However, the limit on g 
necessary to keep the axino mass below the overclosure bound will be somewhat weakened 
with respect to the ordinary case since the trilinear coupling A for g^QiQ2 originates from 
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radiative corrections. 



As a generic situation we consider models which do not fulfill our criterion for a small tree 
level axino mass and possess couplings A = 0(1). The axino mass in such models is given 
by rria = 0(7713/2) (no additional zero modes) or = 0(1713/2/ k) (additional zero modes). 
First of all this mass has to be larger than the mass of the ordinary LSP in the MSSM to 
allow for a decay of axinos. Second, these decays have to occur before the LSP decouples 
which translates into a bound on rria of typically ~ 100 GeV for fpg = 10^^ GeV if the 
decay channel into top and stop is possible. Otherwise an even stronger bound rria ~ TeV 
is required. We see that these generic models are significantly constrained by cosmological 
considerations, however, a final decision depends on details of the model like the exact axino 
mass, the sfermion masses, the PQ-scale etc. Therefore it might be interesting to study models 
which put further constraints on these parameters like e. g. supersymmetric unified models 
incorporating an axion. For a PQ-scale in the upper half of the allowed range fpq > 10^^ 
GeV the lower bounds on rria and the saxion mass (which both increase quadratically 
with fpq) become very stringent and it might be difficult to construct viable models. In this 
context, having axino masses larger than 7713/2 in models with additional zero modes might be 
an interesting option. 
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